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Abstract. The peculiar radial velocities of eight clusters of galaxies, as 
determined from a variety of distance indicators, are consistent with a 
bulk flow not greatly different in direction and magnitude from the veloc- 
• ity vector of the Local Group relative to the 

\ microwave background radiation. The bulk flow vector is directed to 

/ = 266° ± 18°, 6 = 17° ± 15° with a magnitude of 788 ± 113 km s^i. 
N-body simulations of gravitationally induced peculiar motions from a 
set of over 1000 clusters of known redshift give reasonably good agree- 
ment for the peculiar velocities of the Local Group and Virgo both in 
magnitude and direction. Most of the accelerating mass lies within about 
^SJ ■ 40/i~^ Mpc but more distant clusters are necessary to fully account for 

On ' the magnitude and direction. 
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' 1. Introduction 

O ' It has been more than 25 years since the first peculiar velocity of a cluster 

+l3 . of galaxies was reported (Gudehus 1973), more than 20 years since peculiar 

^ I velocities were first referenced to the microwave background radiation (MBR) 

(Gudehus 1978), accurately measured for the first time by Smoot, Gorenstein, 
and Muller (1977), and more than 10 years since Dressier et al. (1987) first 
^ ■ reported that galaxian peculiar velocities were correlated in direction. The first 

H \ part of this time period, however, was marked by the fixed idea that the universal 

expansion was "quiet" (Visvanathan and Sandage 1977 and references therein; 
Tammann 1984). Tammann (1984), for example, claimed that my 1978 peculiar 
velocity for Virgo was "impossible in size and sign" . These early contrary claims 
have been explained (Gudehus 1989b, 1991) as due to i) a neglect to reference 
redshifts to the MBR frame and ii) a, variety of systematic effects. The current 
view is that peculiar velocities exist at some level out to 15,000 km s~^, but with 
some observers finding significantly smaller values at greater distance (Riess et 
al. 1995; Giovanelli et al. 1998) than others (Gudehus 1989b; Lucey el al. 1991; 
Bothun et al. 1992; Mould et al. 1993; Hudson et al. 1999); Willick (1999), 
and one group (Lauer & Postman 1994) finding a significantly different vector 
direction for the most distant clusters. 

In this paper I will use relative distance moduli taken from my previous 
studies, solve for a preferred set of revised cluster peculiar velocities, and com- 
pute an average direction for the cosmic fiow. I will then, from an N-body 
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simulation based on masses assigned to clusters of known redshift, calculate 
predicted peculiar velocities for the LG and Virgo. 

2. The Data 

The data sets used in this study comprise m* magnitudes (Gudehus 1989b) for 
ten clusters, reduced galaxian radius parameters, r^, (Gudehus 1991) for four 
clusters, and relative magnitudes of seventeen selected clusters computed as 
weighted averages from several independent distance indicators (Gudehus 1995) 
(including first-ranked magnitude and richness-corrected nuclear magnitudes). 
The m* and Vg data were combined and fitted by nonlinear least squares to 
the redshift-magnitude diagram with go assumed to be 0.5, luminosity evolution 
assumed to be 0.75 magnitudes per unit z, and all redshifts corrected to the 
Cosmic Background Explorer (COBE) results for the LG's motion relative to 
the MBR (Smoot et al. 1991). The nearest five clusters were given zero weight 
in the fit since their peculiar velocities produce substantial deviations from the 
curve and are nonrandomly distributed. The standard deviation about the curve 
fitted to the remaining points is 0.064 mag. Peculiar velocities were derived from 
the selected clusters in a similar way and include three duplicates from the first 
set, i.e., Virgo, Fornax, and Hydra I, which were given zero weight in the fit. 
Although the standard deviation about the fit is larger, i.e., 0.11 mag, than for 
the first set, the preferred peculiar velocities of the duplicates are taken from the 
second set since in that set each cluster's point represents a weighted average 
of several independent distance indicators (three for Fornax and Hydra I, and 
six for Virgo). Table Q lists for each cluster, the peculiar velocity observations. 
Column 1 gives the identity, column 2 the redshift velocity relative to the MBR, 
column 3 the distance modulus relative to Virgo consistent with the peculiar 
velocities, and column 4 the radial peculiar velocity, respectively. 

Table 1. Cluster Peculiar Velocity Data. 



Cluster 


(km s~"^) 


Am 


(km S~^) 


Virgo 


1473 


0.00 


387 ± 66 


Ursa Major 


1140 


0.25 


-79±111 


Fornax 


1428 


0.26 


202 ± 78 


Hydra I 


4027 


2.44 


665±167 


Perseus 


5179 


3.62 


-541±399 


A1367 


6880 


3.71 


883±200 


Coma 


7185 


4.06 


165±214 


A2199 


9081 


4.71 


-335±297 



Figure |^ shows the radial peculiar velocity vectors projected on the best 
fitting plane passing through the eight clusters. The pole of this plane is in the 
direction a = 6^*7, 6 = 16°, or / = 198°, 6 = 6° (1950.0 equinox). The angle 
between the pole and the MBR apex is 77°. The angle between the pole and 
the supergalactic pole is 31°. 
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Figure 1. Peculiar radial velocity vectors of eight nearby clusters as 
projected on the best-fitting plane through their positions on the sky. 
The weighted average of the vectors, excluding the Local Group, is 
directed towards / = 272°± 18°, h = 48°± 15°, and differs by 19° from 
the direction of the Local Group's motion relative to the MBR. 



A bulk flow vector, solved for by nonlinear least squares, is directed to 
/ = 266° ± 18°, 6 = 17° ± 15° with a magnitude of 788 ± 113 km s'^. The mean 
error of the fit is 136 km s~^, and the vector differs by 15° from the direction 
of the LG's motion relative to the MBR. The average depth is 4549 km s~^ 
and the average weighted depth is 2214 km s~^. Bulk flow tangential velocity 
components added to the observed radial components give the space vectors 
plotted in Figure |2|. In Figure |3| I have plotted on a Hammer- Aitoff projection 
in galactic coordinates, the directions of the bulk flow vectors of this study and 
several others. 



3. Simulations 



In the simulations presented here I consider the gravitational accelerations due 
to clusters of galaxies for which redshifts are presently known. In the first step 
of the simulation calculation, clusters within a given redshift range are selected 
from a catalogue of all known clusters. For each cluster the coordinates and 
distances based on a given Hubble parameter, Hq and deceleration parameter, 
go, are computed. The distances are then reduced to values appropriate for a 
young universe and the model universe allowed to expand. At each uniform time 
step, the acceleration from each cluster on the LG and Virgo are computed, and 
the velocities are computed from the integrals. Ten free parameters are used in 
the simulation. Besides Hq and go, there are the mass of a richness 2 cluster, 
^A2', the factor to be applied to the number of galaxies in a richness cluster in 
describing a cluster of unknown richness, F; the power describing how cluster 
mass scales with richness (see below), P; the age of the universe at the beginning 
of the calculation, ti; the minimum distance cilloweci between clusters, t^wixi] the 
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Figure 2. Peculiar space velocity vectors of eight nearby clusters as 
projected on the best-fitting plane through their positions on the sky. 
The tangential components are taken from a bulk velocity vector fitted 
to the radial components by nonlinear least squares. The bulk flow 
vector is directed towards / = 266°ib 18°, b = 17°ib 15°with a magni- 
tude of 788 ± 113 km s~^, and differs by 15° from the direction of the 
Local Group's motion relative to the MBR. 
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Figure 3. The directions, on a Hammer- Aitoff projection in galactic 

coordinates, of the observed bulk flow velocity vectors of the present 
study and those of Dressier et al. (1987) (D), Hudson et al. (1999) 
(H), Lauer and Postman (1994) (LP), and Wilhck (1999) (W), all with 
respect to the microwave background radiation (MBR), whose direc- 
tion is shown as well. Also shown is the direction of the LG's motion 
as derived from an N-body simulation with 1009 clusters of known 
redshift. 
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number of time steps, N; and the minimum and maximum redshifts used in 
selecting the clusters from the catalogue, z^am and Zmax- 

The cluster mass, Ai2, was taken from Abell (1974) and Kent and Gunn 
(1982) as 2.3 x W^^Mq, and F was arbitrarily set at 0.1. Richer clusters contain 
a proportionally greater amount of light which exceeds the ratio of increase in 
the number of their member galaxies, because the luminosity functions of richer 
clusters extend to brighter magnitudes (Gudehus 1995). The amount of mass 
in richer clusters increases faster than the rate of luminosity increase because 
richer clusters possess higher mass-to-luminosity ratios (Davies, et al. 1980). I 
model cluster mass as a function of number of galaxies of a given richness by 
the relation A4r = Mi^rl'^-if ^ where n is the average number in a richness 
class. An estimate of P can be obtained from the velocity dispersions in clusters 
of two different richnesses. Taking values for the Coma cluster (Kent and Gunn 
1982) and cluster A1689 (Gudehus 1989a), I derive P = 2.1. Based on Larson 
(1969), the age at the beginning of the calculation is taken to be 1.0 x 10^ y. 
The minimum cluster separation is assumed to be 1 Mpc, but this parameter 
rarely comes into play. The Hubble constant is taken to be 70 km s~^Mpc~^. 
A number of time steps equal to 1000 gives sufficient accuracy, and the redshift 
limits are taken as and 0.3. 

A simulation with the above values includes 1009 clusters of which 1 is of 
richness 5, 10 of richness 4, 61 of richness 3, 131 of richness 2, 257 of richness 1, 
358 of richness 0, and 190 of unknown richness. The resultant peculiar velocities 
and galactic coordinates of the vectors are shown in Table |2|. Increasing P to 
2.8 gives better agreement with the observed velocities of the LG and Virgo's 
radial peculiar velocity. The angular deviation of the LG's direction from the 
MBR is 27°. The direction for the simulation with P = 2.8 has been included 
in Figure ^. 

Table 2. Cluster Peculiar Velocity Simulations. 



p 


VhG 
(km s~^) 


^Virgo 

(km s^^) 


^r, Virgo 

(km s^"*^) 


longLG 

longvirgo 


latLG 


2.1 


600 


584 


162 


304° 


23° 










316° 


3° 


2.8 


626 


712 


348 


308° 


33° 










317° 


17° 



Figure |^ shows a smoothed false color plot of the mass assigned to the 
clusters closer than z = 0.07 and which lie within 45° of the best-fitting plane 
through the eight nearby clusters. Since the best fitting P is greater than the 
value based on known cluster masses, i.e., 2.1, increasing amounts of mass are 
contained outside the richer clusters. Figure ^ shows how the model velocity 
for the LG depends on upper redshift limit. Most of the accelerating masses lie 
within about AOh'^ Mpc {z = 0.013), but more distant clusters, such as those 
at 150/i^^ Mpc {z = 0.05) contribute as well. Convergence to a steady value 
occurs only after a distance of about 450/i~^ Mpc {z = 0.15) is reached. 
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Figure 4. A contour plot of mass in the N-body model used for sim- 
ulating peculiar motions, which is associated with clusters of known 
redshift that are closer than z = 0.07, and which lie within 45 degrees 
of the best-fitting plane through the eight nearby clusters. The loca- 
tions, representing 366 clusters, were smoothed with a gaussian of 1232 
km s~^. The Local Group is at the center. The peaks represent extra- 
galactic mass concentrations (EMACS) made up of very rich clusters 
and their associated dark matter. For example, the peak at the left 
at about z = 0.05 includes the richness 4 cluster Abell 3558 and the 
richness 3 cluster Abell 3559. The Perseus cluster is the small peak 
just to the right of center. The MBR flow is to the left. 
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Figure 5. The magnitude of the pecuhar velocity of the Local Group 
as a function of the upper redshift limit of the N-body model. Most 
of the peculiar velocity is accounted for by clusters within about 40 
/i~^Mpc {z = 0.013). An additional contribution at 150 h~^Mpc {z = 
0.05) arises from clusters which include Abell 3558 and Abell 3559. 
Convergence to a steady value is reached after about 450 /i~^Mpc {z = 
0.15). 



4. Summary 

Both the magnitude and direction of the observed average radial peculiar velocity 
are consistent with the set of eight clusters participating in the same flow that 
the LG moves with. Gudchus's (1995) observation that when the measured infall 
to Virgo of 52 km s^^ is used to compute the LG's undisturbed velocity vector, 
the projected radial velocity precisely matches Virgo's measured radial peculiar 
velocity, demonstrates that at least Virgo participates in the flow. Other clusters 
in that study also showed evidence of flow participation. 

Simulations for gravitationally induced peculiar motions in an expanding 
universe with model parameters set equal to reasonable values drawn from the 
literature give acceptable agreement with the observed velocities of the LG and 
Virgo. The sparseness of the cluster catalogue prevents extending the calculation 
to more distant clusters at this time. 
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